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ON THE FORM AND POSITION OF THE SEA-LEVEL AS DEPENDENT ON 
SUPERFICIAL MASSES SYMMETRICALLY DISPOSED WITH RE- 
SPECT TO A RADIUS OF THE EARTH'S SURFACE. 

By Mr. R. S. Woodward, Washington, D. C. 

[continued from vol. ii, page 103.] 

9. Some special values of the integrals (22)-(25), and the corresponding 
potentials (20) and (21), are worth deriving. These values are 

a. For a point of the disturbed surface at the centre of the disturbing mass 
a = o, and (24) gives 



/j = 2 sin i/9j ^n = ^ sin |/9; (26) 



o 



y = \TT 



and the corresponding potential is 

Vi = 4r^hpn sm\^. (27) 

b. For a point at the border of the disturbing mass a = /9, and h.ence froiu. 
(24) or (25) 

/, = /j = 2 arc sin I sin \^ sin y I 

r=o (28) 

a result which is reached quite as readily from (22) or (23). The corresponding 
value of the potential is 

l\= V, = 4r,hp^. (29) 

c. For a point of the disturbed surface 180° from the centre of the disturbing 
mass a = 71, and {25) gives 

\^ 
^ ^ ^ rsin^|^^os^2V^2 
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J I + CDs'* i/3 tan V2 



2^ 



(30) 



■n — 2 cos ^/5 



Ij^ arc tan tan (cos ^fi tan y^) 

o 

= 2;rsin^|j8; 
and the potential is 

Fj = 8 r(,/2/>ff sin' ^ji. (31) 

d?. Suppose ji = JT; in other words, let the attracting mass cover the whole 
sphere ; then (24) gives 

A = ^, (32) 

and we have the following well known approximate value for the potential of a 
spherical shell for a point on its exterior or interior surface, viz : — 

4ro¥7t. (33) 

This result follows also from (27) or (29), if we make /3 = ;r. 

10. In deriving the expressions (20) and(ji) for the potential of the disturb- 
ing mass, it was assumed that a sufficient degree of approximation is attained if 
quantities of the orders «/r„, vjr^, and upwards are neglected. The grounds of 
this assumption need to be examined with some care. For this purpose we shall 
derive the exact expressions, in form at least, for the potential of the disturbing 
mass at its centre, at its border, and at 180° from its centre. A comparison of 
these exact values with the approximate values given in (27), (29), and (31) will 
show the order of approximation of (20) and (21). 

We will first write down the exact expression for the potential of a spherical 
shell of uniform thickness and density, for a point within its bounding surfaces. 
This expression will be u.seful as a check on formulas relating to partial shells. 
Let the radius of the interior surface of the complete shell be r„, the uniform 
thickness k, the uniform density p, and the distance of the attracted point from 
the centre of the shell r^ -\- v; then the potential is * 

r^_ / r I ^ W^ I hv 1 ^ -j , , 

(/-4Vv;r|_n-^^^^_^^^ 2 (r„ -F .;) /^ + 2 (r„ + t/) r, e{r,^v)r,h\ ^34J 

1 1. Let the notation be the same as that used heretofore; i. e. let r be the 

See Price's Calculus, Vol. III. p. 299. 
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radius-vector of any point of the attracting mass, r' that of any point of the dis- 
turbed surface; /> the density, h the thickness, and ^ the angular radial extent of 
the mass ; and 6 the angular distance between the attracted and attracting points, 
the origin being at the centre of the sphere. For points of the disturbed surface 
lying above the undisturbed surface, r will be less than r' = r„ -j- z/ over the range 
r„ to r' , and r will be greater than r' over the range r' to r„ -f- k. Bearing these 
facts in mind, it follows that the exact value of the potential of the disturbing 
mass, for the point where its axis pierces the disturbed surface, is 

sin dd 



V= 2f>7: 



■Pdrf 



o 

h /3 



|/(r^ + r'"^ — 2rr' cos ff) 
sin d M 



r „, C s\n tfdil 

+ ^PV ' ^'' J ,/{^ + r'''-2rr'cost)) 

r„+v o 

= 2pnf[v[^r' sin^ ii? + (r' - rf]- {r' - r)] '^ 
+ 2p7zf\_y[4rr' sin^ i/J + (r - r'f]- (r - r')] ~. 



'-0 + " 

By expansion, integration, and reduction this expression gives, to terms of 
the fir-st order inclusive, 

(35) 
V=Ar hm sin liS f I + lhsm\^-h t^(i-si ni^) ^ ] 

This equation, if we make i? = ?r, agrees with (34), as it should, to terms of 
the second order. It differs from (27) by terms which must be small unless h 
and V are very large. In one of the most important applications discussed by 
the writer, A= 10 000 feet, z" = 3000 feet and ^ = 38°. With these values, 
since r„ is in round numbers 21 000 000 feet, the quantity in the brackets of (35) 
differs from unity by less than I / loooo. If /9 = 60°, which is (as shown in the 
sequel) the angular extent of mass required to produce the maximum elevation 
of the disturbed surface at the centre of the mass, the quantity within the brack- 
ets of (35) exceeds unity by less than 1/8000, using the above value of h and v. 

12. Similarly, if the attracted point be at the border of the attracting mass, 
the exact value of the potential is 

V= 2pj^, drj'y'Kr' — rf + 4rr' sin^ 1/9 J • ^^ 
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+ 2pfp drj^i{r - r'f + ^rr' sin^ |^J . dX 
— P^^{r> — r) p dr — />;rj {r — r') ^ dr, 

''o '-0 + ^ 

_ sin^/?cos^>^ 



Introducing this value, the first of the above integrals relative to X becomes 

>^^cos^A+[j^-^]'cos^r 



(r' + r)f^- 



I — sin^/3sin^ A 



.dk. 



Since the numerator of the element-function in this integral is greater than sin ^ 
cos X, and less than sin 3 cos / -f -r— — cos 8, the value of the integral lies between 

o 
Therefore, if Q^ denote a suitable factor, the exact value of the integral may be 
written thus : — 

Likewise represent the exact value of the second integral relative to X by 



(r+r')^\i + Q., 



'r -}- r' 

Then, to terms of the first order, inclusive, the above expression for the potential 
becomes 

/^ (3 - Q2) ^(i + Q,) 



V=4r„/i/i^ 



I + 



+ 



(36) 



t^(6 l-f Q2) _ 5 /i' - 2hv + 2t^ 

4r„/r " >• 8(r, -ft/)/; ^ 

The first term of this agrees with (29). The factors Q-^ and Q,^ lie between 
o and the elliptic of the first species 

cos /? 



)s I r dX 

^ J 1 (i — sin^^sin^/ 
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When fi = \~, or when the attracting stratum covers a whole hemisphere, Qi = 
Q2 =: o, and (36) agrees, as it should, with the half of (34) to terms of the second 
order; and in this case if h = loooo feet and v = 3CXX) feet, the quantity in the 
parenthesis of (36) exceeds unity by I / 4000. When /9 < ^tz, Q^ and Q^ will ex- 
ceed o, and the degree of approximation of (36) will be somewhat greater than 
in the case just discussed. 

13. For a point of the disturbed surface i8o° from the centre of the disturb- 
ing mass, the exact value of the potential is 

ro + h /? 

IT C 1 1 C sin dd 

V= 2p^j ^drj ^^,.^,n^,,,r,^,y) 
ro + i 

= ifm \\r -\- r' — ^/(r^ + r'^ + 2rr' cos ^)] -, dr. 
'■0 
Expanding, integrating with respect to r, and reducing, there results, to terms 
of the first order inclusive. 



V= 8 r.kpn sin^ i/9 [ i + ^|r+^- ] • (37) 



4 ('-o + 2') . 

The first term of this agrees with (31). Using the above values of r, and k, and 
'd = — 1000 feet {v being here intrinsically negative), the factor in the brackets 
of (37) exceeds unity by I / 2600. 

14. The preceding expressions for the potential of the attracting mass, 
namely, F as defined by equations (20) and (21), are sufficient for mo.st of the 
cases they were designed to meet. They possess the obvious advantage of a 
compact integral form. For some purposes, however, it is desirable to have V 
expressed in a series of spherical harmonics, or Laplace's functions. We may 
thereby arrive at equations (20) and (21), by a process differing from that followed 
in § 7, and establish a useful harmonic development of the elliptic integrals /, and 
/j of equations (22) to (25). 

Expressions fulfilling the present requirements may be derived from equa- 
tion (10) by expanding Z)~Mn a series of ascending powers of -^or . Thus 



equation (9) we have 




I if f r^ f'' V 


. . , when r < r'; 


I I r fr'^ Cr'^ ^ 


. . , when r > r'. 



In these equations J'^, P^, 1\ . . . are Laplace's co-efficients of the zero, first. 
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second, etc. order respectively. They are functions of the angular co-ordi- 
nates f), 0', X, X' only. 

Taking now the centre of the sphere of reference as the origin of co-ordi- 
nates, and supposing the line from which 6 and 6' are reckoned to pass through 
the centre of the attracting mass, we shall have for a mass of uniform thickness 
h, all integrations in (lo) independent. For that part of the disturbed surface 
which lies above the undisturbed surface, r' ^r^-\- v will fall between the ex- 
treme values of r, which are r, and r^,-^ h; and hence 

r <ir' for values between r, and r, -f v, 
r > ^ for values between ^o + z' and r„ + h. 
For that part of the disturbed surface lying below the undisturbed surface, r > r'. 
In both cases the limits of are o and the angular radius /3 of the attracting 
mass ; and the limits of X are o and 2n. Therefore for that part of the disturbed 
surface lying above the undisturbed surface, equation (ro) gives 

v= ,/si„ 6de^dxj[p, + p, [^] + p, [^] V . . . Y^ 

o or, 

+ f,^s\nedejdXj\p, + P, [^] + P, [^'] \ . . . ydr. (38) 

O O r^-lf-v 

Likewise, for that part of the disturbed surface which lies below the undisturbed 
surface, equation (10) gives 

^ 27r r„ -1- .4 

V= /> Jsin ddejdXJ[p, + P, [^] + P, [^] % . . . ^dr. (39) 
o o '•0 

For brevity, let the integrals with respect to r in (38) be denoted as fol- 
follows : — 

■^' -J U' J r' - {i + 3) {r, + z>y + ' ' ^^°^ 

ro + A , 

The ambiguous form assumed by (41) when t = 2 is easily shown to be equal to 

(^-0 + ^)^10^.,^. 
Similarly, let the integrals in (39) with respect to r be denoted thus : — 
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/".• = ^^j^' [''o-'^-^' - (ro + ^)-''-''] . (42) 

Substituting these equivalents in (38) and (39), they become respectively 

V =^ p f Sin dddf[(J„+J'„)P, + i/,+J\)P,+ . . .^dk, (43) 

o o 

(3 2ir 

V = f>f sin dddfW'o P„ + J'\ P,+ . . .^dX. (44) 

o o 

Now it is known from the theory of spherical harmonics * that the general 
value of Pi is 

Pi ^/i (cos d).fi (cos 0') + terms multiplied by cosines of multiples of {k — k'). 
But since. 

27r 

I cos i {X — k') dX = o, 

o 

we have to deal only with the first term oi Pi. If we write for brevity // = cos 6, 
the function of cos 6 or cos d' involved in this first term, is defined thus : — 

Mcose)=fi{f.) = ^-^-^-^--~..-^?^^. 

The following important relation exists between any three consecutive values of 
/i (ji), viz. : — 

/i (/') = ^^^ M-i (/^) — '-^/-2 (/^) • (45) 

The remaining integrals in (43) and (44) are therefore of the form 

/3 2TV I 

/i (cos 0')f/i (cos d) sin ddefdX = 2 k/^ (/i')//- (/^) d/i. 

o o COS 3 

I 

Let F,(;9)=J/,(/i)^/y. (46) 

cos/3 

The known value of this integral is 

J'i (/?) = ^7^ [/i - 1 (cos ^) -/. + , (cos ^)] . (47) 

The values of Fin (43) and (44) may now be written thus, replacing 0' by 

*See Heine, Handbuch der Kugelfunctionen, or Todhunter's Treatise on Laplace's, Lamfs, and 
BesseVs Functions. 
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by a, 0' or a being the angular distance of the attracted point from the centre of 
the atti acting mass: — 

V= 2pn 'l{J, + J',) f, (cos a) F, (§) , (48) 

j = o 
for points of disturbed surface above undisturbed ; 

F= 2pJlJ'\ A (cos a) F, 09) , (49) 

/ = o 
for points of disturbed surface below undisturbed. 

15. The values of Fin the last two equations are exact, but their applica- 
bility to practical problems is limited by the slow convergence of the .series in the 
second members. In most cases with which we are concerned, r„, r^^ v and 
r^-\- h are nearly equal, and hence the convergence depends almost wholly on the 
convergence of the functions of a and ^. 

If we neglect terms of the order ^/r„, f/rj, and upwards, equations (48) and 
(49) become identical and equivalent to (20) and (21). To show these facts, we 
expand the second members of (40), (41), and (42), and obtain, to terms of the 
second order, 

J, + J', = r,h [i —^^[{i -l)h- 2iv + (/ + i) J }] (50) 

(i — ^)h iv~\ 



J"i = r,h^^ 



n> 



(50 



Hence, to terms of the first order, we have * 

Ji+A=J"i = r,h, 
and (48) and (49) become 

V = 2r„ hftnlfi (cos a) Fi (^) . (5 2) 

j=o 

Now if cos <p = cos a cos ^ + sin a sin 6 cos {). — )-') , 

2Tt Zfi (cos a) Fi (d)= \ I -77 J— I — \ • 

,•=0 ' • vf'' J J .^(i _ 2 cos ^ -f l) 

o o 

But this integral, as shown in § 7, is equivalent to /, or /j of (22) or (23) according 
as a is less or greater than /3. This establishes the equivalence of (52) with (20) 

*This inference from (50) and (51) does not appear to be quite satisfactory. For large values of i, 
Ji-\-J'i, andy'j are less than r^h; they are each zero for j'^oo . The influence of the too great factor 
r,, h in the higher terms of (52) is counteracted, however, by the small factors yi (cos o) Fi (/3) in those terms. 
The order of approximation of (52) is evident from its equivalence with (20) and (21), whose order of ap- 
proximation is investigated in W 10-13. 
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and (21), and furnishes the following development in polar harmonics of the 
elliptic integrals in (22) to (25) : — 

o 

-if (S3) 

o 

16. In order to show the forms of the second members of equations (52) and 
(53)1 we give below the first four values of/j (cos «) and Fi{^). The series may 
be easily extended by means of the relations in (45) and (47). 

/o(cosa)= I, 

/i (cos a) = cos a, 

/j (cos a) =: ^ (3 cos'' a — i), 

fg (cos «) = ^ (5 cos* a — 3 cos a) ; 

^o(/5)=i-cos/3, 

^.(^) = i(i-cos=/5), 

F,(/9)-i(cos^-cos='/3), 

^3(/3) = i(6 cos^/9 - 5 cos''/9 - i). 

By means of these values eqnjation (52) may be written thus : — 

2 sin^ i/3 + Y ^°^ "■ s'"' ? 1 

+ J (3 cos" a — I) sin* /3 cos ;9 
+ tV (5 cos'« — 3 cos a) (6 cos" ^9 — 5 cos* /? 



V= 2r^hpit 



L+ . . . J 



17. In case the disturbing mass is as large as the supposed ice mass of the 
glacial epoch, the attraction of the re-arranged free water on itself may be ap- 
preciable. To determine the exact effect of this attraction would be a work of 
great difficulty even if we had the requisite information, namely, an accurate 
knowledge of the complicated shapes of the continents and sea-bottom. But we 
may determine an effect which will in general exceed the probable actual effect 
by supposing the whole surface of the earth covered with a film ot water free to 
assume the proper form for equilibrium under the given forces. To fit this 
ideal case, formula (6) may be modified in the following manner: — 

Let the potential Fin (6) be replaced by V -\- liV, where J F is the potential 
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due to the re-arrangement of the water. Likewise, replace the constant V„ by V^ 
-\- JV^. Then if Av denote the corresponding change in v, equation (6) gives 






(55) 



Now V -{- Jv may be expressed by a series of Laplace's functions (which are in 
this case polar harmonics), thus : — 

v-^Jv = r^c, {Z, + Z, + Z,+ . . .), 

in which c„ is a constant of small numerical value; and hence, denoting the den- 
sity of sea water by />„ , we have, as shown by Laplace*, 



JV= ^r^c.pjz (Z, -f \Zy + \Z^ -f 



2? -f I ' ^ 



). 



V has already been expressed in equation (54) in a series of Laplace's functions. 
Denoting these functions in (54) for brevity by 1^, Fi, Fj, etc., 



V=2r,hpTi(Y,^ Y^^Y^^- . 



). 



Substituting these values of V, AV, and 7/ -)- Av in (55), there results, if we make 
the obviously permissible substitution. 






Ar^Cn-^ 



the following equation : — 



+ Vo^.[l-|g 



[ 



3 P" 



:k+^k)> 



3/.f (n-c/j 

3^ ^ K 



+ . . . 

+ r,c,Z, { I - -4- -^^-"1 -3^-^ 
^^ " « • I 2t + 1 (>„) ^ p.. 

According to the theory of spherical harmonics, we must have in this equation 
the sums of the harmonics of the same order separately equal to zero, which 



*Mecanique Cileste, Book III. Chap. 2. 
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amounts to placing each line in the equation equal to zero. Hence we have 



J = 



-V'i 



{' 



Vi 



(/. 



3 



■1. 



2t + I />,„ 






(56) 



This equation expresses the total effect of the disturbing mass in altering the 
sea-level, v being the effect which would result if the ocean were an infinitely 
rare fluid, and Jv being the increase over v which would result under the assumed 
conditions. Obviously v and Jv may be expressed separately. Thus : — 



Jv • 



3I". 



3^0 ]^ 



[to be continued.] 
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ON THE LAST CLAtJSE OF ART. I4O OF GAUSS S THEORIA MOTUS. 

The clause mentioned is as follows : — 

"Calculo examinando sequentes quoque formulae inservire possunt, quarum 
tamen derivationem non ita diflicilem brevitatis caussa supprimimus : — 

_ a sin Q" — I') _ b sin {I" — I) sin {d' — a) sin { I' — I) 
°~ R R' - sind' ^ W ' 

_ R' sin d' f/cos^os /?" 



R" sin d" • sin {AD' — S) sin s' ' 



ubi f/exprimit quotientem -r 



5 



Tsm (t + r') 



(art. 138. aequ. 10)." 



sin {d' — a) cos {d' — a) 

Notwithstanding the words " non ita difficilem," the derivation of these equa- 
tions is not generally readily seen by the student reading the Thcoria Motus for 
the first time. 

The first equation is obtained from equation III of the same article, by sub- 



